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ABSTRACT 

This paper presents a Novel Multi-Objective Evolutionary Algorithm used for Optimal Reactive Power (VAR) 
transmit. The optimal reactive power transmit difficulty is formulated like a non-linear constrained multi-objective 
optimization difficulty where the real power loss as well as the bus voltage deviations are to be reduced at the same time. 
A novel technique called NSGA-II as well as Strength Pareto Evolutionary Algorithm (SPEA) based approaches are 
planned to handle the difficulty as an exact multi objective optimization difficulty. A hierarchical clustering algorithm is 
compulsory in SPEA to provide the decision maker to maintain controllable pareto -optimal set. Furthermore the results are 
compared with NSGA-II and found SPEA is slight superior than NSGA-II. 

The objective difficulties are formulated by using GA along with different constraints as well as studied utilizing 
projected techniques. Each and every constraint was satisfied as well as the results obtained for IEEE 30-bus systems are 
presented, illustrating the effectiveness of the projected approach. 

KEYWORDS: Reactive Power and Devices, Constraints, MOEA, SPEA, NSGA-II 

INTRODUCTION 
Evolutionary Algorithms 

Evolutionary Algorithms (EAs) such as evolution strategies and Genetic Algorithms have turn into the technique 
of preference for optimization difficulties that are too multifarious to be solved utilizing deterministic methods such as 
linear programming or gradient (Jacobian) techniques. EAs need a little information about the difficulty being solved, and 
also they are simple to put into practice, strong, and also inherently parallel. To determine a certain optimization difficulty, 
it is sufficient to need that one is capable to calculate the objective purpose for a specified set of input parameters. 
Because of their universality, simplicity of operation, and also condition for parallel computing, EAs frequently take fewer 
time to finish the optimal resolution than gradient techniques. One of the goals of an ideal multi -objective optimization 
process is to determine as many pareto -optimal resolutions as possible. Since an EA works with a population of 
resolutions, a population of pareto -optimal resolutions can be captured in one single simulation run of an EA. If we can 
attain this it will remove recurring utilize of a single objective optimization technique for determining one different 
pareto -optimal resolution in each run. This will also remove the requirement of any such parameters like weight vectors, 
e vectors, target vectors and different others. 
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Some of the other advantages of having evolutionary algorithms are that they need extremely small information 
about the difficulty being answered, fewer vulnerable to the shape or continuity of the pareto front, simple to put into 
practice, strong and could be implemented in parallel environment. 

The most important motivation for utilizing EAs to resolve multi -objective optimization difficulties is because 
EAs deal at the same time with a set of possible solutions which permits us to determine some members of the 
pareto-optimal set in a single run of the algorithm, instead of having to carry out a series of separate runs as in the case of 
the conventional mathematical programming methods. In addition, EAs are fewer vulnerable to the figure or continuity of 
the pareto front, whereas these two issues are a real concern for mathematical programming methods. 

• Non-Dominated Sorting GA: The technique is recommended by Goldberg which is based on some layers of 
categorizations of the individuals. Before selection is performed, the population is ranked on the basis of 
non-domination. To continue the diversity of the population, these classified individuals are mutual with their 
replica fitness standards. Then this cluster of classified individuals is ignored and another layer of non-dominated 
individuals is considered. The procedure keeps on until all the individuals in the population are classified. 
A multi-objective optimization difficulty has a number of objective functions which are to be reduced or 
maximized at the same time. In this study a recent method called "Strength Pareto Evolutionary Algorithm" is 
utilized to reduce the objective functions. 

• Non-Dominated Sorting GA-II: This is the modified edition of the NSGA, called NSGA-II, which is further 
resourceful, utilizes superiority and also a crowed comparison operator that continues diversity without specifying 
any additional parameters. The NSGA-II doesn't utilize an external memory as in the previous algorithms. Instead 
its superior mechanism consists of combining the most excellent parameters with the most excellent off spring 
attained. 

• Strength Pareto Evolutionary Algorithm: Zitzler and Thiele in 1998 projected a best evolutionary algorithm, 
which they called the Strength Pareto Evolutionary Algorithm (SPEA). This algorithm starts with a randomly 
generated population 'P 0 ' of size 'N' and an unfilled external population 'P t ' with highest capability on 'N t \ 
In any generation, the most excellent non-dominated explanations are copied to the external population set 6 P t \ 
Subsequently the dominated results in the personalized external population are found and removed from the 
external population. In this manner, previously found elites which are currently dominated by latest elite results 
obtain removed from the external population. What remains are only the most excellent elite results. If this 
process is continued over several generations, there is a hazard of the external population being overcrowded 
along with non-dominated results. In order to control this overcrowding of population above 'N t ', grouping 
methods is utilized. 

Once the latest elites are protected for the subsequent generation, the algorithm then turns to the current 
population and also utilizes genetic operators to determine a latest population. In this the primary step is to assign a fitness 
to each solutions in the population, so that a fitness 'Si' to each member i of the external population first. The strength Si is 
proportional to the number Ni of present population members that an external solution i dominates; 



Si = ni/(N+1) 



eq. (a) 



where N is the number of population in the current population. 
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Thereafter the fitness of a current population member j is assigned as one more than the sum of the strength values 



The addition of one makes the fitness of any present population member T 0 ' to be more than the fitness of any 
external population member P t . this technique of fitness assignment suggests that a solution with a smaller fitness is better. 

SPEA- Algorithm 

Step 1: Find the most excellent non-dominated set of T 0 ' and copy these solutions to P t , where 
P 0 , is initial population and P t is the external pareto set with N and N t respectively. 

Step 2: Find the most excellent non-dominated solutions of P t of the personalized population P t and remove all 
dominated solutions. 

Step 3: If IP t l> Nt, utilize a grouping method to decrease the size to N t . Otherwise keep P tj unchanged. The 
resulting population is the external population P t +1 of the subsequent generation. 

Step 4: Allocate fitness to every best solution i G P t+1 by utilizing equation, eq (a) then allocate fitness to every 
population member j G P 0 by utilizing equation eq (b) 

Step 5: Apply a binary tournament selection with these fitness values in a Minimization sense, a cross over and 
mutation operator to generate the latest population P 0+ i of size N from the mutual population P t+1 +P 0 of size (N t + N). 

Step 3 and step 5 results in the latest external and present population which are then processed in the subsequent 
generation. This algorithm continues until stopping criterion is met. 

NON-DOMINATED SORTING GENETIC ALGORITHM - II (NSGA-II) 

NSGA (Non-Dominated Sorting in Genetic Algorithms) is a popular non -domination based genetic algorithm for 
multi-objective optimization. It is a very effective algorithm but has been generally criticized for its computational 
comp- lexity, lack of elitism and for choosing the optimal parameter value for sharing parameter a share. A customized 
version, NSGA-II was developed, which has a improved sorting algorithm, incorporates elitism and no sharing parameter 
requires to be chosen a priority. 

The modified version of the NSGA, called NSGA-II, which is more resourceful, utilizes elitism and a crowed 
comparison operator that keeps diversity without specifying any additional parameters. The NSGA-II doesn't utilize an 
external memory as in the previous algorithms. Instead its elitist mechanism consists of combining the best parameters 
with the best off spring obtained. 

The primary step of this technique is to sort the population P according to non -domination. Different algorithms 
can be utilized for this task. It is important to reiterate that any two members from the similar class cannot be said to be 
superior to one another with respect to all objectives. The total number of classes, represented as 'p' in the following 
equation, depends on the population P and the under lying problem. 



Once the classification job is over it is clear that all solution in the first set that is, all i G Pi, belong to the most 




of all external population which weakly dominate j; Fj = 1+ iePt ^J 



eq (b) 



P = U"i=iPj 
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excellent non-dominated set in the population. The second most excellent solutions in the population are those that belong 
to the second set or all members of P 2 and so on. Obviously the worst solutions are those belonging to the final set or all 
members of P p , where p is the number of different non-dominated set in the population. 

The fitness assignment process starts from the first non-dominated set and consecutively proceeds to dominated 
sets. Any solution i of the first non-dominated set is allocated fitness equal to Fi = N. This exact value of N is utilized for a 
particular function. Because all solutions in the first non-dominated set are uniformly significant in terms their proximity to 
the pareto optimal front relative to the present population. 

In an NSGA, the sharing function technique is utilized to supply the diversity of the solution. That is, for each 
solution i in the front P i9 the normalized Euclidian distance dy from another solution j in the same front is calculated as 
follows. 



It is significant to note that the sharing distance is calculated with the decision variables, instead of the objective 
function values as in an MOGA. Once these distances are calculated, they are utilized to calculate a sharing function value 
utilizing the subsequent equation with a = 2. 



The sharing function obtains a value between 0 and 1, depending on the distance dy. Any solution i who has a 
distance better than a share from the i th solution contributes nothing to the sharing function value. After all IP II sharing 
function values are considered, they are added together to compute the niche count nci of the i th solution. The niche count, 
in some sense represents the number of solutions in the neighbourhood of the i th solution. If there is no other solution exists 
within a radius of a share from a solution, the niche count for that solution would be one. 

The ultimate task is to decrease the fitness of the i th solution by its niche count and get the shared fitness value. 
F/ = Fi / nci. This method of degrading fitness of a solution which is crowded by several solutions helps importance the 
solutions residing in less crowded region. This process completes the fitness assignment process of all solutions in the first 
front. In order to continue the second front, we note the minimum shared fitness in the first front and then assign a fitness 
value slightly lesser than this minimum shared fitness value. 

Crowding Distance 

Once the non-dominated sort is concluded the crowding distance is assigned. Since the individuals are selected 
based on rank and crowding distance all the individuals in the population are assigned a crowding distance value. 
Crowding distance is assigned front wise and comparing the crowding distance between two individuals in different front 
is meaningless. The crowding distance is considered as follows: 

• For each front F h n is the number of individuals. 

• Initialize the distance to zero for all the individuals i.e. Frfdj ) = 0, 
Where j corresponds to the j th individual in front F t . For each objective function m 
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o Arrange the individuals in front F t based on objective m i.e. I = sort(F im ). 

o Assign infinite distance to boundary values for each individual in F ( . i.e., I(di) = oo and I(d n ) = oo 

o For k = 2 to (n - 1) 

I(d k )=I(d k )+[I(k+l).m-I(k-l).m]/f m max -f m min 



Where I(k).m is the value of the m th objective function of the k th individual in i. 
Selection 

Once the individuals are arranged based on non -domination and with crowding distance allocated, the selection is 
carried out utilizing a crowded comparison operator (<n). The comparison is carried out as follows based on: 

• Non-domination rank prank i.e. individuals in front Fi will have their rank as p rank = i. 

• Crowding distance F t (dj) p <n q if - p rank < q rank 

• Or if p and q belong to the same front F t then i*Xd p ) > Fj(d q ) i.e. the crowding distance should be more. 
The individuals are selected by utilizing a binary tournament selection with crowd comparison operator. 

Genetic Operators 

Real-coded GA's uses Simulated Binary Crossover (SBX) operator for crossover and polynomial mutation. 
Simulated Binary Crossover 

Simulated binary crossover simulates the binary crossover observed in nature and is give as below. 

c ltk =|[(i-A)Pi,*+(i+A)P2,*] c 2Jt =|[(i+A)A,*+(i-A)P2,*] 

where C i k is the i th child with k th component, p ik is the selected parent and fi k (> 0) is a sample from a random 
number generated having the density, 

p [P) = 0.5(7 +¥\ P ( " ] =^M)-^, 

Wc )H ' if 0 < p < 1 P ifp>l 



This distribution can be obtained from a uniformly sampled random number u between (0, 1). r| c is the distribution 
index for crossover. That is 

i /?(«) = ! r- 



J3{u) = (2w) (7+,) [2(1 - 

Polynomial Mutation 

c k =p k +(pt ~p[)S k 

Where c k is the child and p k is the parent with Pk being the upper bound on the parent component, * k is the 
lower bound and * is small variation which is calculated from a polynomial distribution by using 
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S k -<2r.)*-* -1, if rk<05 S k =l-[2(\-r k W" + \ ifrk > 0 . 5 

Where r k is a uniformly sampled random number between (0, 1) and n " > is mutation distribution index. 

Recombination and Selection 

The offspring population is now combined with the current generation population and selection is performed to 
set the individuals of the next generation. Since all the previous and current best individuals are added in the population, 
elitism is ensured. Population is now sorted based on non -domination. The new generation is filled by each front 
subsequently until the population size exceeds the current population size. If by adding all the individuals in front Fj the 
population exceeds N then individuals in front Fj are selected based on their crowding distance in the descending order 
until the population size is N. And hence the process repeats to generate the subsequent generations. 

PROBLEM FORMULATION 

The optimal VAR dispatch problem is to optimize the steady state performance of a power system in terms of one 
or more objective functions while satisfying several equality and inequality constraints. Generally the problem can be 
formulated as follows. 

Objective Functions 

Real Power Loss (Pl): In this case sum of the losses of all branches in the system is to be minimized. This results 
in the reduction of active power generation in the system. This in turn saves both the generation cost and creates a higher 
generation reserves. This objective is to minimize the real power loss in transmission lines that can be expressed as; 

•A = Pi = E& W? + vj - ivy, .- 8j ,)] 

k=\ 

Where nl is the number of transmission lines, g k is the conductance of the k th line, and 7 ' are the 

voltages at the end buses i and j of the k th line respectively. 

Voltage Deviation (V D ) 

This objective is to minimize the deviation in voltage magnitudes at load buses that can be expressed as; 

V ref 

Where NL is the number of load buses; k is the pre-specified reference value of the voltage magnitude at the 

yref 

k th load bus. A is usually set to be 1.0 pu. 
Problem Constraints 

• Equality Constraints 

These constraints represent typical load flow equations as follow 
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MB 

/= 1,2 NB 

\ j.> 

> =1 f = 1,2 NB 

Where NB is the number of buses, P G and <2 G are the generator real and reactive power respectively, P D and Q D 
are the load real and reactive power respectively, Gy and #y are the transfer conductance and susceptance between bus i and 
j respectively. 

• Inequality Constraints 

• Generation constraints 

Generator voltage V G and reactive power outputs Q G are restricted by their lower and upper limits as; 

vr<v r <vr 

U/ 0/ 0/ , i=l,2, NG where NG number of generators 

er^G^GT ,_,, 2 , ng 

• Transformer constraints: 
Transformer tap T settings are bounded as; 

ymin < < 7^ max 

1 ' , i = 1,2, NT where NT is the no. of transformers. 

• Switchable VAR sources constraints: 

The Switchable VAR compensations Q c are restricted by their limits as; 
O min < O < n max 

&d *£ci Ma ? i = \ 9 2 9 , . .NC, where NC is the number of capacitors. 

• Security constraints: 

These includes the constraints of voltage at load buses V L as follows; u ~ u ~ u , i = 1,2, NL 

where NL is no of Load buses Aggregating the objectives and constraints, the problem can be mathematically formulated 
as a nonlinear constrained multi -objective optimization problem as follows; 

Minimize [P L (x,u), VD(x,u)] subjected to; g (x, u) = 0 h (x, u) < 0 

Where, x is the vector of dependent variable consisting of load bus voltages V L , and generator reactive power 
outputs Q G . There are a set of control variables, which can be modified to satisfy the load flow equations. Hence X can be 
expressed as; 

u- is the vector of control variables consisting of generator voltages V G , transformer tap settings T, and shunt 
VAR compensations Q c . Hence u can be expressed as ; u ~ ^ V ° x " Vq ^ ' ^ ' " TwT 9 "' ^ ^ 



www.tjprc.org 



editor@tjprc.org 



66 



G Kalidas Babu, Rajesh Kumar Samala & K Swathi 



g— is the equality constraints, 
h— is the inequality constraints. 

A general multi-objective optimization problem consists of a number of objectives to be optimized simultaneously 
and is associated with a number of equality and inequality constraints. It can be formulated as; Minimize f t (x), 1=1, 
2, N oh] 

Subject to: gj (x) = 0 j = 1 M, constraints and h k (x) <= 0 

k = 1 K, constraints 

Where £ is the i th objective functions, x is a decision vector that represents a solution, and N oh j is the number of 
objectives. For a multi-objective optimization problem, any two solutions x 1 and x 2 can have one of two possibilities, one 
covers or dominates the others or none dominates the other. In a minimization problem, without loss of generality, solution 

x 1 dominates x 2 if the following conditions are satisfied: V/ e I 1 ' 2 ' N otj } • ft( x ) ^ ) 

3/e{l,2, N ob ,}zf J (x l )<f J (x 2 ) 

If any of the above condition is violated, the solution xl does not dominate the solution x2. If solution xl 
dominates the solution x2, xl is called the non-dominated solution. The solutions that are non-dominated within the entire 
search space are denoted as pareto -optimal and constitute the pareto-optimal set or pareto-optimal front. 

MAIN ALGORITHM (SPEA-METHOD) 

Read all system data including bus data, line data, maximum and minimum limits of Generator bus voltages, Load 
Voltages, Tap Values, Shunt Values, the bus no. at which taps are connected power, etc 

Form Y-bus with present system data 

Generate Initial Population by considering required number of Strings and Population size. 
Decode the strings and evaluate (finding fitness) 
Calculate the actual values of the parameters 

Check the limits of all calculated values and fix to their violated limits 
Run load flow with modified system parameters 
Check all the constraints 

Calculate Objective function ie. P-Loss and Volt deviation as per equation 
Find Nondominated solution from the above Multiobjective function 
Initialize External Pareto Set to Null 
Copy these Nondominated solution to External Pareto Set, 

Check Nondominated Solutions in the External Pareto set and Delete all dominated Solutions from the Pareto set 
Check the size of the Pareto set, if it is greater than pre-defined size, reduce the size by means of Clustering. 
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• Calculate the Fitness Values of individuals in both Pareto set and the population as follows; 

• Assign fitness 'Si' to each member i of the external population first. The strength Si is proportion to the number 
(ni) of current population that an external population T dominates, S A = ni / (N+l) 

• Thereafter the fitness of a current population member 'j' is assigned as on more than the sum of the strength 

Z st 

values of all external population members which weekly dominate j, Fj = 1+ l&Pt * l -J 

• Arrange the chromosomes in descending order of their fitness. 

• Apply genetic operators 

• Uupdate the chromosome matrix. 

• Check the convergence conditions, if satisfied go to step 20 else to step no 4. 

• Stop and Print the result. 

• Stop and print the data 

MAIN ALGORITHM (NSGA-II) 

In this case also run the Load Flow and get all the necessary data i.e. real power loss and voltage deviation. Take 
these as the population for the algorithm. In this, real coded GA is employed. The initialized population is sorted based on 
non-domination. The fast sort algorithm is described as below for each 

• For each individual p in main population P do the following 

• Initialize Sp = 0. This set would contain all the individuals that is being dominated by p. 

• Initialize np = 0. This would be the number of individuals that dominate p. 

• For each individual q in P 

if p dominated q then add q to the set Sp i.e. Sp = Sp U {q} 
else if q dominates p then increment the domination counter 

• for p i.e. np = np + 1 

• if np = 0 i.e. no individuals dominate p then p belongs to the first front, Set rank of individual p to one i.e p rank = 
Update the first front set by adding p to front one i.e F t = F t U {p} 

• This is carried out for all the individuals in main population P. 

• Initialize the front counter to one. i = 1 . 

• Following is carried out while the i th front is nonempty i.e. Fi ^ 0 

• Q = 0. The set for storing the individuals for (i + l) th front. 

o for each individual p in front Fi 
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o for each individual q in Sp (Sp is the set of individuals dominated by p) 
o nq = nq-1, decrement the domination count for individual q. 

o if nq = 0 then none of the individuals in the subsequent fronts would dominate q.Hence set qrank = i + 1. 
Update the set Q with individual q i.e. Q = Q U q. 

• Increment the front counter by one. 

• Now the set Q is the next front and hence Fi = Q. 

o Now the crowing distance is calculated as below 

• For each front Fi, n is the number of individuals. 

• Initialize the distance to zero for all the individuals i.e. Fi(dj ) = 0, where j corresponds to the j th individual in 



• For each objective function m 

o Sort the individuals in front F t based on objective m i.e. I = sort(Fi,m). 

o Assign infinite distance to boundary values for each individual in Fi. i.e. I(dl) = oo and I(dn) = oo 
o fork = 2to(n-l) 



o The selection is carried out using a crowded comparison operator (<n). The comparison is carried out as 
below based on 

o Non-domination rank prank i.e. individuals in front Fi will have their rank as p rank = i. 
o Crowding distance Fi(dj) 

• p<nqif - Prank < qrank 

• Or if p and q belong to the same front Fi then Fi(dp) > Fz(dq) i.e. the crowing distance should be 
more.The individuals are selected by using a binary tournament selection with crowed-comparison-operator. 

o Apply Polynomial Mutation (refer above theory) 

o Recombination and Selection (refer above theory) The offspring population is now combined with the 
current generation population and selection is performed to set the individuals of the next generation. 

• Repeat the step till it converge. 

• Stop and Print the result. 
CASE STUDY 

In this study, the proposed method was tested on the standard IEEE 30-bus, 6-generator test system to investigate 
its effectiveness. The single line diagram of the IEEE 30-bus system is given in appendix- 1 figure 14. The system has six 
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generators at bus 1, 2, 5, 8, 11, and 13 and four transformers with off -nominal tap ratio in lines 6-9, 6-10, 4-12, and 27-28. 
The lower and upper level voltage magnitude limits are 0.95 and l.lpu at all generator buses 2, 5, 8, 11, and 13 
respectively and 1.05pu for the remaining buses including slack bus 1. The lower and upper limits of the transformer 
tapings are 0.9 and l.lpu respectively. In this case study, Fast decoupled load flow was run as preliminary step. From the 
load flow we get the system conditions and losses in the system. Also we can find the voltage deviations in the load buses 
from the reference value. In this case study the reference value of voltage is taken as 1 pu. This real power loss and voltage 
deviation is taken as the input to the algorithm. First run the load flow with base case and then take the readings (Pl oss and 
VD) for optimization using linear combination of P Loss and VD. 

IEEE 30 BUS SYSTEM LINE DIAGRAM 




Figure 1 

IEEE 30-Bus System: No. of buses: 30 No. of lines: 41 
CASE-1 (SPEA): Parameters Selected for Study 

No of buses=30 No of lines = 41 Generator buses =6 Tap changing transformers =4 No. of shunts =2 
Chromosome Size=250 Total String Length=108 String length for P=07 (each) String Length for V=07 (each) 
String Length for Transformer=4 (each) String Length for shunts=.04 (each) 
RESULTS 

Base Case Voltage Deviation (pu) (VD alone as the objective function) 



-VO_TAGE VARJAT! C T-4 — 



1 6 13 20 



Figure 2: Voltage Deviation Vs Generation 

Base Case Power Loss Variation (pu) (PL alone as the objective function) 
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Figure 7: Power Loss Variation Vs Generation Linear Combination of 
Real Power Loss and Voltage Deviation 30-Bus System 
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Figure 3: Power Loss Vs Voltage Deviation (250-Population, 150- Generation, 35- Pareto) 
Table 1: Result from the Linear Combination of P-Loss and VD 



Pair* lit* ters 


Ee.^r PL 




Powier L-oss 


0.05419 


0.0577 


Voltage Deviation 


0.87S3 


0.2577 


VG-1 


1.07 


1.05 


VG2 


1.0424 


1.01 


VG-5 


1.04 


1.034 


VG-S 


1.02S7 


1.0143 


VG-11 


1.0(532 


1.045 


VG-13 


1.064 


1.0523 


Tap-1 


1.0956 


1.0564 


Tap -2 


0.9872 


1.00 


Tap- 3 


1.08 


0.976 


Tap-4 


1.002 


0.9671 


Shunt 


Cl=0.01, 


Cl=0.04, 




C 2=0. 05 


C2=0.02 



No. of Non-dominated Solutions in Single Run case (Linear combination) 
POWER LOSS (pu) 
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Figure 4: Convergence of Power Loss in Single Run Case 

Secondly the problem has been applied to the proposed approach (SPEA) with the above mentioned parameters 
(i.e. PLoss and VD). Where objective functions, power loss and Voltage deviation on the Load buses are considered 
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simultaneously for the optimization purpose. The problem was run with SPEA and the diversity of the pareto -optimal set 
over the trade-off surface was found. The result obtained by the proposed method is shown in table 2 for different 
conditions. The result is also compared with linear combination method and showing that, the proposed method have better 
diversity characteristics and well distributed over the entire trade-off surface. 

Pareto Optimal front in 30 Run case (SPEA) 

Table 2: Results Shows the Best Values for 250-pop, 150 Generation 



Pri lTiiiieteis 


Brat FL 


Best Vcl 


Power Loss 


0.051724 


0.05529 


V olta ge De viatio n 


0.88143 


0.242 


VG-1 


1.05 


1.035 


VG2 


1.044 


1.0386 


VG-5 


1.0354 


1.0254 


VG-S 


1.03456 


0.9938 


VG-1 1 


1.0646 


1.0268 


VG-1 3 


1.09256 


1.0445 


Tap-1 


1.0998 


1.055 


Tap -2 


0.95 


0.8823 


Tap -3 


0.9625 


0.9754 


Tap -4 


0.95 


0.9243 


Shunt 


C 1=0. 04, 


C 1=0. 03, 




C2=0.05 


C2=0.01 



Voltage Deviation (Pu) 

*** PARETO OPTIMAL FRONT*** 
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Figure 10: Power Loss Vs Voltage Deviation (250-Population, 150- Generation, 35-Pareto) 



Convergence of Power Loss convergence of Voltage Deviation 

Power Loss Variation in Single Run Case Volta § e Aviation in Single Run Case 




Figure 5 

CASE-2 (NSGA-II) Parameters Selected for Study 

No of buses =30 No of lines =.41Generator buses =6Tap changing transformers =4 No. of shunts =2 

Chromosome Size=200 Generation=500 String length for V=06 (real coded) String length for Transformer = 
04(each), 4*4 String Length for shunts=04 (each), 2*4 
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In this method simulation was done with unadjusted fast decoupled load flow to get the power loss and voltage 
deviation. Take this objective function as input to the proposed method and run NSGA-II algorithm for optimal dispatch of 
reactive power. 

*** NONDOMt NATE D SOLUTION "* 
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Figure 6: Non Dominated Solutions NSGA-II Voltage 

Deviation (pu) Voltage Deviation Vs Power Loss 



Table 3: Results Shows the Best Values for 200-Pop, 500 Generation 
(The Best Result from 30 Separate Runs Case) 



FVirri lit- ter?s 


Brat PL 


BestVd 


Power Loss 


0.05209 


0. 05S5 


Voltage Deviation 


0.7S2 


0.2021 


VG-1 


1.06 


1.05 


VG2 


1.042 


1. 1 


VCt-5 


1.03 


1.041 


VG-S 


1.027 


1.026 


VG-11 


1.072 


1.05 


VG-13 


1.0S4 


1.062 


Tap-1 


1.09 


1.05 


Tap-2 


0.9572 


1.00 


Tap- 3 


1.09 


0. 9S 


Tap-4 


1.01 


0. 97 


Shunt 


Cl= 0.03, 


Cl=0.04, 




C 2=0. 05 


C2=0.01 



Table: 4 Comparisons of Results Obtained by Different Methods 







Comb ma lion 








Best PL 


Best va 


Best PL 


Best Vd 


Best PL 


Best vd 


Best PL 


Best Vd 


Po^ver Loss 


0.0559S 


0.0872 


0.05419 


0.0577 


0.051724 


0.05529 


0.05209 


0.0585 


YD 


0.889 


0.301 


0.8783 


0.2577 


0.88143 


0.242 


0.782 


0.2021 


VG-1 


1.08 


1.04 


1.07 


1.05 


1.05 


1.035 


1.06 


1.05 


VG2 


1.076 


1.04 


1.0424 


1.01 


1.044 


1.0386 


1.042 


1.1 


VG-5 


1.052 


1.038 


1.04 


1.034 


1.0354 


1.0254 


1.03 


1.041 


VG-S 


1.043 


1.03 


1.0287 


1.0143 


1.03456 


0.993 S 


1.027 


1.026 


VG-11 


1.071 


1.042 


1.0632 


1.045 


1.0646 


1.0268 


1.072 


1.05 


VG-13 


1.061 


1.052 


1.064 


1.0523 


1.09256 


1.0445 


1.084 


1.062 


Tap-1 


1.074 


1.062 


1.0956 


1.0564 


1.099S 


1.055 


1.09 


1.05 


Tap-2 


1.045 


1.05 


0.9872 


1.00 


0.95 


0.9823 


0.9572 


1.00 


Tap-3 


0.9821 


1.025 


1.08 


0.9^6 


0.9625 


0.9^54 


1.09 


0.98 


Tap-4 


1.0 


0.987 


1.002 


0.96^1 


0.95 


0.9243 


1.01 


0.97 


Shunt 


Cl=0.04 
C2=0.03 


CI =0.03 
C 2 =0.0 5 


CI =0.01 
C 2=0. 05 


CI =0.04 
C2=0.02 


C 1=0.04 
C2=0.05 


CI =0.03 
C2=0.01 


CI =0.03 
C 2=0. 05 


Cl=0.04 
C2=0.01 



• Base Case Voltage Deviation (VD alone as the objective function) 



Impact Factor (JCC): 5.9638 



Index Copernicus Value (ICV): 3.0 
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GEE N EE RATI ON S 



Figure: 6 Voltage Deviation Vs Generation 

CONCLUSIONS 

• In this paper, a Genetic Algorithm based approach has been presented and applied to multi -objective VAR 
dispatch optimization problem. The problem has been formulated as multi-objective optimization problem with 
competing real power loss and bus Voltage deviation. 

• A new technique called NSGA-II and SPEA methods are used the study. In SPEA a hierarchical clustering 
technique is also applied to keep manageable Pareto -optimal set without destroying the characteristics of the 
trade-off front. 

• The result shows that the proposed methods are efficient for solving multi -objective VAR dispatch problems 
where multiple solutions can be found in one simulation run. Since the proposed approach does not impose any 
limitations on the number of objective functions, its extension to include more objective is a straightforward 
process. 

• Also this method is useful for nonlinear problems with more objective functions. The result shows that the SPEA 
is little efficient than NSGA-II for solving multi-objective VAR dispatch problems. But in terms of simulation 
time, real coded NSGA-II is far better than SPEA. 
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